Recursive calculation of the microcanonical density of states 
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For a classical system of noninteracting particles we establish recursive integral equations for the 
density of states on the microcanonical ensemble. The recursion can be either on the number of 
particles or on the dimension of the system. The solution of the integral equations is particularly 
simple when the single-particle density of states in one dimension follows a power law. Otherwise it 
can be obtained using a Laplace transform method. Since the Laplace transform of the microcanon- 
ical density of states is the canonical partition function, it factorizes for a system of noninteracting 
particles and the solution of the problem is straightforward. The results are illustrated on several 
classical examples. 
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I. INTRODUCTION 

In statistical mechanics, an isolated macroscopic sys- 
tem at equilibrium is characterized by its microcanonical 
entropy, given by the Boltzmann formula 



S — kjilnD, . 



(1) 



where fee is the Boltzmann constant and fl the number 
of accessible microstates in phase space, i.e., the number 
of equiprobable microstates compatible with the fixed 
values of the macroscopic external parameters such as 
the volume V, the number of particles N, the energy E, 
etc. For a classical system in D dimensions, the phase 
space has 2 s dimensions corresponding to the compo- 
nents of the generalized momenta p and coordinates q 
such that s — DN. Each microstate corresponds to a 2s- 
dimensional cell in phase space with a volume /i", where 
h is the Planck constant. The accessible volume is such 
that 



E <H{p,q) <E + SE, 



(2) 



where the arbitrary uncertainty 6E <C E affecting the 
value of the total energy is introduced in order to obtain 
a 2s dimensional accessible volume, thus leading to a 
nonvanishing value of il. 

The number of accessible microstates is proportional 
to SE and can be written as 



Q = nNiE)SE, 



(3) 



where the coefficient of proportionality nN{E) is the mi- 
crocanonical density of states which plays a central role 
in the statistical physics of isolated systems. 

The calculation of the density of states amounts to de- 
termine the volume enclosed by the hypersurface of con- 
stant energy E = H{p, q) in the 2s dimensional phase 
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space, a quite difficult task in general. In textbooks the 
properties of the microcanonical ensemble are often illus- 
trated by the study of some classical systems involving 
noninteracting particles for which the Hamiltonian is sep- 
arable and the constant energy surface simple enough. In 
the case of the ideal gas one needs the area of a hyper- 
sphere and for a collection of harmonic oscillators the 
area of a hyperellipsoid which can be transformed to a 
hypersphere through rescaling. The area of the hyper- 
sphere is usually obtained by comparing the values, cal- 
culated either in Cartesian or in spherical coordinates, of 
the integral over the 2s-dimensional space of exp(— r^), 
where r is a radius vector (see, e.g., Refs ^M)- Other 
independent particle problems have also been treated^ii 
by making use of Dirichlet's integral formulai^ 

The purpose of the present paper is to propose an- 
other approach to noninteracting particle systems, based 
on a recursion relation for the microcanonical density of 
states. Assuming that the single-particle density of states 
is known, the particles are added one by one. Summing 
their contribution to the total density of states, a recur- 
sive integral equation is obtained which can be solved 
using a Laplace transform method. When the single- 
particle Hamiltonian is itself separable, involving a sum 
of contributions from each of the D dimensions, the same 
recipe can be used to evaluate the single-particle density 
of states in D dimensions. 

The recursive approach to the calculation of the N- 
particle microcanonical density of states is presented in 
Sec. II. Some illustrative examples are discussed in Sec. 
III. Summary and conclusion are given in Sec. IV. 



II. MICROCANONICAL DENSITY OF STATES 

We consider an isolated classical I?-dimensional sys- 
tem containing N noninteracting distinguishable parti- 
cles. Let Np{E) denote the p-particle integrated den- 
sity of states, i.e., the number of microstates with en- 
ergy smaller than _E > in a system with p parti- 
cles. The p-particle density of states, np{E), is such that 



dNp{E) = np{E)dE. 

The number of microstates with energy smaller than E 
in a system with p-\-l particles is obtained by first com- 
bining the microstates with energy smaller than E — e'vn 
a system with p particles with the microstates on the en- 
ergy interval [e, e + de] associated with the added particle 
and then integrating over e. Doing so, one obtains 



Np+i{E)= f Np{E-e)ni{ 
Jo 



e)de. 



(4) 



Taking the derivative of both sides and noticing that 
A'j,(0) = 0, leads to a recursion relation for the j^particle 
density of states 



^p+l 



(E) 



dN.p+i{E) 
dE 



np{E - e)ni{e)de , (5) 



The expression on the right is a convolution product thus, 
taking the Laplace transform of both sides, one obtains 



Zp+l{s) 



ip+i{E)dE 



^-sE 



dE / np{E - e)ni{t)de 



ni{e)e~ 



Zp{s) zi{s). 



s{E- 



'>np{E-e)dE 
(6) 



Iterating this result, the Laplace transform of the N- 
particle density of states is given by 



zn{s) = [zi{s) 



N 



(7) 



when the particles are identical. Otherwise the power is 
replaced by a product over the different contributions. 
The A^-particle density of states is obtained by taking 
the inverse Laplace transform 



UNiE) 



1 

27ri 



■y-'rioc 



e^ ZN{s)ds. 



(8) 



where 7 is a vertical contour in the complex plane chosen 
so that all singularities of zn{s) are on the left of it. 

One may notice that with 1/s = k^T, zjf{s) is the 
A^-particle canonical partition function and Eq. ([T]) sim- 
ply states that for noninteracting particles the iV-particle 
canonical partition function factorizes and is given by the 
A'^th power of the single-particle canonical partition func- 
tion zi(s). 

Let us assume that the single-particle density of states 
varies as a power of the energy 



ni{E) = CiE'''-\ E>0. 
The Laplace transform is then given by 



(9) 



zi(s) = Ci / e 
/o 

gK,l 



— sE pKi — 1 



dE 



e-*t«i-i^+- 



dt = Ci 



r(«;i) 



(10) 



where r(a::) is the Euler Gamma function. For the A^- 
particle system we have 



zn{s) 



[CiTX^)] 

„Nki 



N 



(11) 



The inverse Laplace transform can be deduced from 
Eq. (fTO)) and leads to the TV-particle density of states 



njv(-B) = [CiT{Ki) 



N 



r(iVKi) ■ 



(12) 



In this simple case one may also proceed as follows: 
Using Eq. ([9]) and the power-law expression of the p- 
particle density of states 



np(E) = CpE^"-^ 
in the integral equation ([Sj leads to 



(13) 



np+,{E) = CiCp f e^'-'iE^er^-'de 
Jo 

= CiCp£;'"i+'"''^i f t''^-'^{i-t)'^''-'^dt 

Jo 



(J fj r(Ki)r(Kp) ^^_^j^^^_i 

^ T{k.i + Kp) 



(14) 



The integral in the second line is the Euler Beta function^ 
which is written in terms of Gamma functions in the last 
line. Comparing with the form of Up^i (E) resulting from 
Eq. (J13p we deduce the following recursion relations for 
the amplitude and the exponent: 



Cp+i — CiC, 



T{Ki)r{Kp) 



Kp+l = Kl+ Kp. (15) 



^r(Kl -I- Kp) 

Iterating these relations, one easily obtains 

[C,T{k,)]p 

C„ = —Z7-, ^ , Kp = pKl , (16) 



yp 



r{pKi) 



and the form of the iV-particle density of states given in 
Eq. (fT2l) is recovered. 

For a mixture with Na particles of type a and Ni, par- 
ticles of type 6, we have a straightforward generalization 
of Eq. (fT^ , namely 



i^Na,NbiE) = / nN^{e)nN^iE ~ e)de 
Jo 

= [CariK^r^iariK.r^-— — — r-^ . (17) 

In order to calculate the microcanonical density of 
states, one needs the form of the single-particle density 
of states. 

When the single-particle Hamiltonian is the sum of 
equivalent (kinetic and potential) contributions coming 
from each of the D dimensions of the system, the num- 
ber of dimensions plays the same role as the number of 



particles did previously and a recursion relation similar 
to Eq. ([5]) can be written for the single-particle density 
of states in q + 1 dimensions 



,(«+!) 



{E) 



n[^\E - e)7i^^\e) de (18) 



where the upper index now refers to the dimension of 
the system while the lower one refers to the number of 
particles. 

As before the convolution product leads to a form sim- 
ilar to O 



SD) 



z^is) = [zi^^(.) 



(1)/ 



1-D 



(19) 



for the Laplace transform of the single-particle density of 
states in D dimensions. 

When the single-particle density of states in one di- 
mension behaves as 

n(i)(^)^cf'£;«i"-\ E>Q, (20) 

the single-particle density of states in D dimensions reads 

D £;^'^5"-i 



SD) 



'C\E)^ C[''T{k\'>) 



.(1)^ 



T\Dn'^^^ 



(21) 



in analogy with Eq. (112]) . 

When the system is anisotropic with Di dimensions 
[i ~ a, b) for which the one-dimensional single-particle 
density of states is Ui (E) with Laplace transform zj (s), 



(22) 



Eq. (|T9|) has to be replaced by 

4^-^^)(.) = [z(")(.)]^»[z('')(.)]^^ 



III. SOME EXAMPLES 

A. Boltzmann Ideal gas 

We consider an ideal gas of N point particles of mass 
771, confined inside a cubic box of volume Vjj = L^ . The 
single-particle Hamiltonian is 



Hip) 



P 
2m 



(23) 



In the two-dimensional phase space of the single-particle 
problem there are two segments of constant energy, E, 
such that p{E) = ±y/2mE and < q < L. Hence, as 
shown in Fig. [l] the microstates with energy lower than 
E are located inside a rectangle with sides 2\/2mE and 
L. Since each microstates corresponds to an area h, the 
Planck constant, we obtain 



n[^\e) = ^V2^, n[^\E) 



L 2m 



E ' 



(24) 



(2mE) 



P 

1/2 







-(2mE) 



1/2 



L 



FIG. 1: Rectangular domain of accessible microstates with 
energy lower than E in the phase space of a free particle with 
mass m, confined on a segment with length L. 



for the one-dimensional single-particle densities of states. 
Thus Cf ^ == V2^L/h, k[^'^ = 1/2 and Eq. ^ imme- 
diately gives the single-particle density of states in D 
dimensions: 



.(D) 



(E) = Vd 



h' J r{D/2) 



(25) 



The microcanonical density of states of the ideal gas 
with N particles of mass to in a volume Vd follows from 
Eq. dni) with Ki = D/2, Cir(Ki) = VD{2TTm/h^)D/^ , 
according to Eq. (PS)) . so that 



'■N 



(E) = V, 



N 



f2T:m\'""^ E^D/^-^ 



^ \ h^ J r{ND/2) 

and the integrated density of states is given by 



^/2^to\^^/" g^^/2 
>N y^J- "D y 1^2 ) T{ND/2 + l) 



AD) 



N'->iE) = V, 



(26) 



(27) 



in agreement with the result obtained by standard 
methods.— 1^1^ For undistinguishable particles the densi- 
ties in Eqs ^ and ^ must be divided by iV!, the 
number of permutations between the N particles. 

With an ideal gas consisting of Na particles of mass 777,0 
and Nb particles of mass 7T7f,, such that Na + Ni, = N, 
Eq. ini) leads to 



''Na.Nt 



(E) - V, 



N 



ND/2 



Na Nt 

m„ rrii 



D/2 eND/2^1 



r{ND/2) 

(28) 

for the microcanonical density of states of the ideal gas 

mixture in D dimensions. The integrated density of 



states is given by 



ND/2 






N, 



D/2 E^^/'^ 



T{ND/2 + iy 

n ('^^ 

in agreement with Ref. 16|. A division of the densities by 
NalNbl is needed when the two species consist of undis- 
tinguishable particles. 



B. Classical harmonic oscillators 




FIG. 2: Elliptic Domain of accessible microstates with energy 
lower than E in the phase space of a one-dimensional har- 
monic oscillator with mass m and angular frequency uj. The 
semiaxes are a = \/2mE and /3 — y2E/rn/Lo. 

We consider a collection of N noninteracting classical 
harmonic oscillators in D dimensions which for simplicity 
are assumed to be isotropic. 

For a single oscillator in one dimension the Hamilto- 
nian reads 



H{p,q) 



P 
2m 



moj 



(30) 



The line of constant energy E in the two-dimensional 
phase space is an ellipse with semiaxes a — \2/mE and 
/3 = ^j2Elra/uj enclosing an area 7ra/3 = 2ttE/uj shown 
in Fig. [2] Each microstate corresponding to an area h, 
we obtain the integrated denstity of states as 



n[^\e) 



27r 



-E, 



,(1) 



{E) 



27r 
hiij 



(31) 



.(1) 



Ji) 



Thus C{' = 2TT/huj and k\' = 1. According to Eq. ^ 
the density of states for one oscillator in D dimensions is 



AD) 



'-)h£ 



E 



D-l 



{D^l)\ 



(32) 



so that Ki = _D, Cir(Ki) = {27r/huj)^, and, according 
to Eq. p^ . the microcanonical density of states for N 
oscillators in D dimensions reads 



'AT 



(E) 



2n\ 



ND 



E 



ND-l 



{ND~ 1)! 



and one obtains 



N^n\e) 



2^\ND ^ND 



huj 



iND)l 



(33) 



(34) 



for the integrated density of states in agreement with the 
known result i^iiSiiii The generalization to a collection of 
anisotropic oscillators is immediate; It amounts to re- 
place uj^ by rii^i'^^- 



C. Ideal gas in a gravitational field 

We consider a classical ideal gas with N point parti- 
cles of mass m in D dimensions. The gravitational field 
with acceleration g acts downwards along the vertical di- 
rection. The gas is confined inside a vessel with height I 
in the vertical direction and section So = L^~^ in the 
transverse directions. This is an example for which the 
single-particle density of states is not a simple power of 
E so that we shall use the Laplace transform method. 




FIG. 3: Parabolic domain of accessible microstates with en- 
ergy lower than E in the phase space of a particle with mass 
m, moving in the one-dimensional half-space g > 0, in a gravi- 
tational field with uniform acceleration g directed downwards. 
When the motion is restricted to < q < I, the domain of 
accessible microstates is restricted to the part of the shaded 
region with q < I < E/mg when E > mgl. 

Let us calculate the single-particle density of states 
starting from Eq. (j22p since the system is anistropic. 

For each of the Da = D — I transverse dimensions, 
the single-particle density of states is given by Eq. ([M)) 



so that Ci — \/2mL/h, k^ = 1/2 and according to so that 

Eq. mr 



z^^\s) 



LV2^r(l/2) _ L^/27 



,1/2 



/lsl/2 



(35) 



o /O-^ /'inin( ,1) 

.W...-2V2™/ (™..)^/___^^ 



7v}'''(£;) 



h 

2x/2to 







In the remaining dimension (Df, = 1) we have to cal- 
culate the one-dimensional density of states for a particle 
in a gravitational field. The Hamiltonian reads 



3hmg 



Vtdt 
E^/^-{E - mglf^QiE-mgl)], (37) 



hmg Jmax(_E-mgi,0) 

4\/2mr 



2m 



Q<q<l. 



(36) 



where Q{x) is the Heaviside step function such that 
Q(x) = when x < and Q{x) = 1 when x > 0. The 

The line of constant energy in phase space is the seg- derivative with respect to E leads to 

mcnt of parabola with equation p — zLy/2'my/E — mgq 

in the region < q < I (see Fig. [3]). Thus the number 

of microstates with energy lower than E is obtained by 

dividing by h the surface between the parabola and the p 

axis (restricted to the region where q < I when E > mgl) 

I 



,('') 



riE)- 



2V2m 



hmg 



E^''^-{E-mglY''^Q{E-mgl)\ . (38) 



Taking the Laplace transform of this density of states, one obtains 



z'/\s) 



2V2m 
hmg 



2\/2^r(3/2) 
hmg s^/'^ 



/ E^/^e-'^dE-er'"'ai (E ~ mgl)^/^e-''^^-"'3i)dE 

Jo Jmal 



(1 - e-"™9') 



I mgl 

^/2^ (1 - e-^^fl 



hmg 



„3/2 



(39) 



The Laplace transform of the iV-particle density of states in D dimensions then follows from Eqs ([7]) and p2|) together 
with Eqs dSl]) and ^ 



4^-^.l)(.) = [,W(,)]^(^-l)[,f)(,)]A'^P^ 



Sd 



mg 



N 



2^m\^^/'(l-e'^™sO 



/l2 



N 



„N(D+2)/2 



C \ ^ /o \ ND/2 N /AT\ 



sknigl 



The inverse Laplace transform is obtained using the transformation formula 

for each term in the sum. It leads to the following expression for the density of states 



SD-IS) 
'•AT 



[Ey 



mg 



N 



„ X ND/2 N 

2nm \ x— 



h^ 



b-i)m; 



fc=0 



N\ {E - kmgl) 



N{D+2)/2-l 



T[N{D + 2)/2] 



Q{E - kmgl) , 



(40) 



(41) 



(42) 



and the integrated density of states is given by 

^ ^-- - ^j ^^/2 ^ ^ ^^ (N\ {E - kmgl)^(^+^y^ 

fc=0 

in agreement with the result of Ref. Il2l . 



N: 



{D-l,l)fj^, 



N 



^ ' \mgj \h^ ) ^^ ^ \k) r[iV(D + 2)/2 + l] ^ ^^ 



(43) 



I 

When I -^ oo, only the first term contributes to the can be obtained more directly since the deviation of the 
sum and one recovers the result of Ref. 7. This result single-particle density of states from a simple power law 



coming from the last term, in Eq. (j38l) disappears in this 
hmit. 

Here too, one has to divide the densities by iV! when 
undistinguishable particles are considered. 



IV. SUMMARY AND CONCLUSION 

We have shown how, using the integral equation ([5]), 
one can construct the A^-particle density of states itera- 
tively, starting from the single-particle density of states 
and, in the same way, using the integral equation ([TSl) . 
how to deduce the single-particle density of states in D 
dimensions from the single-particle density of states in 
one dimension for a microcanonical system of noninter- 
acting particles with a separable Hamiltonian. 

The solutions of the integral equations are easily ob- 
tained when the single-particle density of states in one 



dimension follows a power law. This has been illustrated 
on two classical examples: The Boltzmann gas and a col- 
lection of classical harmonic oscillators. In other cases, 
the density of states may be obtained using a Laplace 
transform method. Since the Laplace transform of the 
microcanonical density of states is given by the canonical 
partition function of the system, it factorizes for nonin- 
teracting particles. As a consequence it is easily deduced 
from the form of the single-particle densities of states for 
each dimension as shown for an ideal gas in a vessel of 
finite height under the influence of the gravitational field. 
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